Andreev scattering in nanoscopic junctions at high magnetic fields. 
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We report on the measurement of multiple Andreev resonances at atomic size point contacts between 
two superconducting nanostructures of Pb under magnetic fields higher than the bulk critical field, 
where superconductivity is restricted to a mesoscopic region near the contact. The small number of 
conduction channels in this type of contacts permits a quantitative comparison with theory through 
the whole field range. We discuss in detail the physical properties of our structure, in which the 
normal bulk electrodes induce a proximity effect into the mesoscopic superconducting part. 



PACS numbers: 61.16.Ch, 62.20.Fe, 73.40.Cg 

It is well know that it is possible to fabricate atomic 
size contacs between metallic electrodes in a controlled 
way by means of the mechanically controllable break 
junction technique or the scanning tunneling microscope 
(STM). Indeed, by repeatedly indenting the tip into 
the sample one can achieve a stationary state in which 
a connecting neck between the electrodes is formed. [||^] 
This neck elongates and contracts during the repeated 
indentation following a well defined pattern of elastic 
and plastic steps, which has been neatly measured in 
a combined STM-AFM experiment where conductance 
and forces could be recorded simultaneously || . 

The properties of a given neck can be probed by mea- 
suring the current-voltage characteristic within the same 
experiment, so that the STM serves at the same time 
as a fabrication tool and as an experimental probe of a 
very singular atomic size nanostructure. ^,^] A reason- 
able knowledge of the geometry of the neck which can 
be varied in a well controlled way, is obtained through 
a simultaneous measurement of the conductance during 
the fabrication process Q . The final form of these struc- 
tures, which are successfully fabricated [|[§J^] is a long 
connecting neck jointed on its ends to the bulk electrodes 
whose radius decreases in a smooth way towards a central 
constriction, which can be of atomic size. 

In this experiment, the control of the morphology ex- 
tends over two lengths scales: first, the overall form of the 
neck can be varied at mesoscopic length scales (hundreds 
or thousands of A) by the repeated indentation process, 
and second, the smallest cross section can be varied at 
atomic scales (tens of A) by doing small voltage varia- 
tions on the z- piezotube. Recently new possibilities of 
atomic size contacts have lead to progress on the under- 
standing of some phenomena occurring at a nanoscopic 
level. It has been shown during the last years that lead 
(Pb) is a good material to create this kind of small di- 
mensions systems having the additional advantage of be- 
ing a superconductor below T c = 7.16K. jj| Indeed, the 
transport of current between two weakly linked supercon- 
ductors brings noteworthy information about the contact 



through e.g. the Josephson current or through the mul- 
tiple Andreev reflection mechanism. || In the case of 
a single atom link between two electrodes the authors 
of Ref. H proposed that the effect of the multiple An- 
dreev resonances on the I-V characteristics is a measure 
of the number and transparency of the conduction chan- 
nels through a single atom JToj] . 

In this work we focus on the magnetic field dependence 
of the I-V characteristics of single atom point contacts. 
Indeed, it is well known that superconductors of reduced 
dimensions such as thin films or granular samples remain 
superconducting well above He. |8) As the magnetic field 
penetration depth of lead is about 390A for a bulky sam- 
ple, it is feasible to build connecting necks with smaller 
lateral dimensions with the repeated indentation proce- 
dure. We find indeed that sufficiently long and nar- 
row necks show superconducting features up to fields as 
large as 20 times the bulk critical field of Pb (which is 
0.05T at 1.5K). The structure of subgap resonances due 
to multiple Andreev reflections remains under field and 
our analysis shows in detail how the pair breaking effect 
of the magnetic field, together with the N-S proximity 
effect from the bulk electrodes smears the subgap reso- 
nances. 

We use a stable STM setup with a tip and a sample of 
the same material (Pb) which is brought from the tunnel- 
ing into the contact regime by cutting the feedback loop. 
The I-V curves were taken at 1.5K using a standard four 
wire technique. Great care was taken to shield electri- 
cally the whole setup as RF noise is known to smear the 
subgap resonances in small contacts. The experiment is 
done by gently changing the smallest cross section of the 
neck to make a large number of atomic size contacts at 
each magnetic field without varying the overall form of 
the neck. Indeed, while our setup is sufficiently stable to 
maintain the same neck over a complete magnetic field 
sweep, we cannot maintain the morphology of the neck on 
the atomic level over a large field variation. Nevertheless, 
we could perform small field sweeps of several hundred 
Gauss with a given atomic arrangement and we find the 
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same result, so that the measurement procedure does not 
change the results presented here. The maximal elonga- 
tion of the piezotube, which is 1600 A, limits the overall 
length of the necks. Here we discuss one typical case of 
a neck having a critical field of about 20 times the bulk 
critical field of lead with a magnetic field applied always 
parallel to the long axis. 
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FIG. 1. I-V characteristics of atomic sized contacts in 
the presence of a magnetic field taken at 1.5K. V is normal- 
ized to A at each field (from the top to the bottom: Ao, 
0.89A , 0.78A , 0.7A with the measured A = 1.35meV), 
and / is also normalized to A and to a, the conductance 
above 2A (as expressed in units of the quantum of conduc- 
tance ao = 2e /h). Not all experimental points are plotted 
in order to show more clearly the fits (straight lines) using 
the model explained in the text. The parameters in the lower 
right corner are the transmissions through the different chan- 
nels used to fit the experimental data (see text). The sum 
over the three parameters gives the total conductance a in 
units of o"o- Each line of numbers corresponds to one curve, 
from top to bottom. T is the pair breaking parameter (also 
defined in the text). The upper critical field of bulk Pb is 0.05 
T. 

Fig.l shows a representative choice of measured I-V 
curves of several last contacts before breaking at zero 
field and under field. Each I-V curve, is different at each 



contact and can be well fitted at zero field by the con- 
duction channels model of Ref. ||] (straight lines, up- 
per figure in Fig.l). Accordingly, the experimental I-V 
curves in a last contact show a large variety of behav- 
iors which is slightly changed by varying the morphology 
of the contact at atomic length scales. This is compati- 
ble with measurements that record the conductance at a 
fixed voltage in a large number of contacts and show how 
the conductance of Pb last contacts shows steps compara- 
ble to the quantum unit of conductance, but the average 
value has large fluctuations [jfl). The model of Ref. |i| 
uses (and verifies) the theoretical predictions that the I- 
V curve between two superconductors which are weakly 
linked through a small number of conduction channels is 
highly non linear and varies strongly depending on the 
transparency T of the junction (0 < T < 1; tunnel to 
contact regimes). It turns out, that only one conduction 
channel with a given T is not sufficient to fit the I-V 
curves shown in Fig.l, but that it is necessary to add a 
number N of theoretical curves, each one with a given T n 
between and 1. This was related to the number N and 
transparency T n of channels in each single atom contact, 
where N and the average values of the T n 's depend on 
the element studied. In the case of Pb, this gives N = 3 
with T\ rather opened (most frequently close to 1) and 
T2 j 3 more closed (smaller than 1). In the uppermost part 
of Fig.l the numbers show the experimentally measured 
T n 's for each contact. We will not go into more details 
about this model which is extensively discussed in Refs. 
|,0|l||l|]. In the following, we discuss how to explain 
the data under field. 

We first analyze the influence of the magnetic field in- 
troducing the pair breaking effect in the standard proce- 
dur e fej , as formulated in a wavefunction representation 
]l3| |l~6f It was shown in Q that pair breaking effects 
can be incorporated by modifying the Andreev reflection 
amplitude, a(ui) — u(lu) — y/u 2 (u>) — 1, where u(u>) satis- 
fies Th: 



h r 1 

t = u 1 - r - — 



(i) 



where T = 1/(Atj,(,), T p b is the pair breaking time and A 
is the self consistent superconducting gap including the 
pair breaking effects. This expression is generally valid, 
irrespective of the origin of the pair breaking mechanism 
p"7j . The value of T used in the fittings was assumed to 
be the same for all channels and all I-V curves at a given 
applied field. 

The straight lines in Fig.l shows the fittings which are 
as good as the ones obtained at zero field, provided that 
the pair breaking parameter is introduced. The number 
N and the characteristic values of the parameters for the 
transparency of each channel T n does not vary up to the 
largest fields. T is determined with a precision of about 
20%. 

The values of T explain the magnetic field dependence 
of the gap in the tunneling regime, which we have mea- 
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sured by breaking completely the contact. We therefore 
could follow and fit precisely the predicted influence of 
pair breaking in the structures associated with multiple 
Andreev reflections. This was not possible to do pre- 
viously, as other realizations of multiple Andreev reflec- 
tions (e.g. large point contacts, tunnel junctions with 
microbridges S) involve experimental setups which are 
much more complex than a single atom contact with a 
small number of conduction channels and cannot be mod- 
elled precisely. 

We can gain more insight in the physics of this system if 
we consider the pair breaking parameter in a uniform su- 
perconducting cylinder in a magnetic field which is given 
by [0: 

h e 2 DR 2 H 2 

1~ = MS (2) 

where R is the radius of the cylinder, D = vpl/3 is the 
diffusion coefficient, I is the mean free path, and H is the 
applied field. Following this model, in order to explain 
the values of r used in fig.jlj we need a cylinder of a 
radius which is rather large (R « 450A, taking £ w R, 
note that a smaller value of £ leads to even larger values 
of R) as compared to the usual width estimations for the 
neck presented here or other necks fabricated with the 
same method Q. Clearly, a model based on a simple 
cylinder does not explain the observed behavior, we need 
to take into account that the radius varies as a function 
of z. 

We analyze in the following the order parameter, den- 
sity of states and pair breaking parameter in a neck of a 
varying radius. Indeed, a better agreement is obtained 
if we consider that at a given field, the superconduct- 
ing region is in good contact with the part of the neck 
with larger radius which already became normal, so that 
pair breaking effects arise from the proximity effect of 
this normal region. Assuming that the electronic mean 
free path is smaller than the coherence length, we can 
describe the superconducting properties by the Usadel 
equations H-OT • We parametrize the Green's functions 
in terms of an angle parameter, 0(r,E), where E is the 
energy measured from the chemical potential plj. Set- 
ting % = 1, they can be written as: 

y V 2 + iEsin(9) + |A| cos(0) 

- 2e 2 D\A\ 2 coa(e) sin(0) = (3) 

where A = (Hru$)/2 is the vector potential. We neglect 
the influence of other spin flip and inelastic processes. 
At the boundary of the contact, we have |V0|_r = 0, 
and R(z) determines the geometry of the neck (we ne- 
glect any radial dependencies, and take R(z) < £). We 
also assume that the magnetic field is unscreened within 
the neck. Then, A can be replaced by its average, 

(A 2 )(z) = H . Within this approximation, the vec- 
tor potential enters in eq.|| as giving rise to an effective, 



position dependent, pair breaking time. If we apply eq. 
|3j to a uniform wire, this pair breaking time reduces to 
that in eq.[|). 
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FIG. 2. Superconducting order parameter as func- 
tion of distance, for different applied fields. The con- 
tact region is located at z/£ — (from top to bottom, 
H = 0.1, 0.2, 0.6, 1.0, 1.2, 1.3 T). The magnetic field is applied 
along the long axis of the cones. 



Figure shows the superconducting order parameter, 
for different fields, as function of the position for a typical 
neck modelled by two truncated cones of L = 800 A 
length attached to bulk electrodes, with an opening angle 
of a = 35°. We also take L = 3£, so that £ w 260l. g| 
There is a smooth transition to the normal state as the 
radius of the neck increases. This is further illustrated in 
figure H|, where the density of states is shown at different 
positions for H — 0.2T. For this field, the influence of the 
normal region is felt throughout the entire neck. Even at 
the central region, the gap is significantly rounded. This 
is also observed in the calculated density of states at the 
center shown in Fig. 4. 
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FIG. 3. Electronic density of states for different positions 
of the neck. The field is H = 0.2T. 
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FIG. 4. Density of states at the contact region for different 
applied fields. 

/From the solution of the Usadel equations we can infer 
the amplitude for Andreev reflections at the contact sur- 
face, which is given by i tan[6*(z = 0, E)/2] . This function 
is slightly different from the standard expression used to 
incorporate pair breaking effects in a point contact (see 
equ. Jl| and Ref. |l4|] ) . We have checked that there are no 
appreciable differences in the quality of the fits to the ex- 
perimental data, shown in fig. Q , with a reasonable value 
for f fa 260 - 30fjA. 

In conclusion, we have measured, and analyzed, the 
multiple Andreev scattering resonances of atomic sized 
Pb contacts in the presence of a magnetic field greater 
than the bulk critical field. In this regime, superconduc- 
tivity is restricted to a small neck of mesoscopic dimen- 
sions. We are able to build and control in situ with our 
STM structures which are a unique example of weak links 
of dimensions variable from atomic to mesoscopic length 
scales, opening a new field of studies in nanophysics. We 
present a quantitative comparison of experiment and the- 
ory of pair breaking effects on multiple Andreev reso- 
nances. 
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